THE  APPLICATION  OF  LINEAR  SUPERPOSITION 
TO  BLEICH'S  TWO-BEAM  THEORY 


Jose  Elvio  de  Assis  Rodrigues 


■ 


THE  APPLICATION  CI  LINEAR  SUPERPOSITION 

TO  BLEICH'S  TWO-BEAM  THEOR? 

by 

Jose  Elvio  cle  Assis  Rodrigues 


B.  S.,  United  States  Coast  Guard  Academy 
(1966) 


si ittee  in  p:  itial  iuliillment 

01  THE  REQUIREMENTS  J?OR  ''  '"   EOREES  01 
!rASTER  01  SCIENCE  IN 
NAVAL  ARCHITECTURE  AND  MARINE  ENGINEERING 
AND  MASTEJ?  01  SCIENCE  IN 


_  0  . 


THE  APPLICATION  01  LINEAR  SUPERPOSITION  TO 
BLEICH'S  TWO-BEAM  THEORY 
by 
Jose  Elvio  de  Assis  Rodrigues 


Submitted  to  the  Department  of  Ocean  Engineering  on  May  12,  1972, 
in  partial  fulfillment  for  the  degree  of  Master  of  Science  in 
Naval  Architecture  and  Marine  Engineering  and  Master  of  Science 
in  Mechanical  Engineering. 


ABSTRACT 


The  principle  of  linear  superposition  is  applied  to  Bleich's 
two  general  solutions  pertaining  to  hull-deckhouse  interaction.   The 
resulting  longitudinal  stress  distributions  for  three  separate  mathemati- 
cal models  are  compared  with  stress  distribution  results  obtained 
through  the  use  of  finite  element  techniques.   In  reference  to  Bleich's 
theory,  procedures  are  derived  for  the  determination  of  the  stress  dis- 
tribution at  locations  away  from  the  center  of  the  hull-deckhouse  struc- 
ture.  Analytic  and  finite  element  procedures  are  also  described  for 
the  determination  of  the  deck  stiffness. 

The  principle  of  linear  superposition  is  found  to  be  valid 
for  Bleich's  theory.  Because  the  theory  neglects  shear  lag  effects, 
its  application  is  recommended  only  for  the  center  portion  of  struc- 
tures with  relatively  long  deckhouses. 
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INTRODUCTION 

In  1953  H.  11.  Bleich  presented  a  paper  in  the  "Journal  of  Applied 
Mechanics"  entitled  "Non-Linear  Distribution  of  Bending  Stresses  Due 
to  Distortion  of  the  Cross  Section."   In  this  paper  he  derived  a  via- 
ble analytical  solution  to  the  problem  of  hull-deckhouse  interaction. 
However,  because  of  the  difficulty  involved  in  evaluating  the  two- 
coupled  differential  equations  for  real  life  situations  Bleich  was 
only  able  to  present  two  simplified  solutions.   One  solution  to  the 
differential  equations  was  simplified  by  considering  only  bending 
moment  loads  and  neglecting  the  distributed  loads.   The  other  solution 
considered  only  distributed  loads  and  neglected  bending  moment  loads. 
Although  these  simplifying  assumptions  are  not  justified  in  any  real 
ship,  the  methods  of  solution  are,  in  actuality,  straightforward  and 
rather  simple  to  obtain.   The  proposal  is  that  the  principal  of  linear 
superposition  can  be  applied  to  Bleich's  solution  in  order  to  deter- 
mine the  total  solution.   Even  though  the  title  of  Bleich's  article 
implies  non-linearity,  the  calculated  stress  distributions  found  were 
actually  linear  in  form  with  the  only  discontinuity  or  break  anpearing 
at  the  main  deck  level.   In  order  to  confirm  or  disprove  this  hypothe- 
sis, finite  element  methods  using  the  ICES  STRUDL  program  developed 
at  M.I.T.  will  be  employed  to  find  the  total  bending  stress  distribu- 
tion.  A  comparison  of  the  linear  superposition  of  the  two  solutions 
using  Bleich's  method  and  that  of  the  total  ICES  STRUDL  solution  should 
give  some  indication  of  the  applicability  of  linear  superposition  to 
Bleich's  method.   If  superposition  if  applicable,  a  rather  simple 
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analytical  solution  could  then  be  available  for  the  determination  of 
the  actual  total  bending  stress  distribution. 


CHAPTER  I  -  BLEICH'S  TWO-BEAM  THEORY 

Consider  the  problem  of  two  separate  beans  forced  to  act  together 
by  horizontal  shear  forces  and  vertical  forces  acting  at  the  junction 
of  hull  and  deckhouse.   In  essence,  the  system  consists  of  a  beam 
elastically  supported  by  another  beam  with  a  shear  connection  to  enforce 
equal  strains  at  the  deck  level,  Pigure  1.1  and  1.2.   The  vertical 
forces  are  due  to  elastic  resistance  of  the  deck  framing  or  bulkheads 
against  the  motion  of  the  deckhouse  with  respect  to  the  hul 1 .   Navier's 
hypothesis  for  a  structure  acting  as  a  beam  is  applicable  only  if  the 
cross  sections  of  the  beam  do  not  distort.   Since  the  side  plating  of 
a  deckhouse  is  not  coplanar  with  the  ships  side, the  vertical  deflections 
of  the  two  are  not  necessarily  the  same.   Navier's  hypothesis  that 
plane  sections  remain  plane  and  the  conventional  theory  of  bending  of 
beams  are  not  applicable  for  the  entire  structure,   for  this  analysis, 
however,  Bleich  assumed  that  the  symmetrical  cross  section  consists  of 
two  non-distorting  portions  which  can  move  relative  to  each  other  and 
that  Navier's  hypothesis  is  applicable  for  the  hull  and  deckhouse  each 
by  itself. 

In  the  structure  shown  in  Pigure  1.1,  the  lower  hollow  box  beam 
represents  the  hull  and  is  of  length  L  while  the  upper  box,  the  deck- 
house, is  shorter  and  is  of  length  1.   Both  boxes  are  assumed  to  be  of 
constant  cross  section.   The  cross-sectional  area  and  the  moment  of 
inertia  of  the  deckhouse  and  of  the  hull  are  A1 ,  I.,  and  A„ ,  I«, 
respectively,  and  the  distances  of  the  respective  centers  of  gravity 
from  each  other  and  from  the  deck  are  a,  a  a,  and  3-a. 
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1«   Simplified  Analysis  of  Two- Co 11  Structure 

In  this  analysis  the  assumption  is  made  that  the  deck  A-B,  Mgure 
1.1,  and  its  supports  have  no  stiffness  and,  hence,  will  not  resist 
any  relative  vertical  movements  between  hull  and  deckhouse.   This 
simplification  is  not  justified  for  any  real  ship  system  and  is  dis- 
carded in  the  final  analysis;  however,  some  important  results  are 
derived. 

At  a  distance  z  in  the  free  body  diagram  of  figure  1.2,  the  moment 
and  direct  forces  in  the  deckhouse  and  hull  are  M  ,  N-  ,  and  M_ ,  N_ , 
respectively,  with  positive  moments  producing  compression  at  the  top 
of  the  respective  units.   Direct  forces  N.  and  N~  are  positive  if  they 
create  tension.   The  external  loads  acting  to  the  left  of  the  section 
have  a  moment  M.   A  shear  force  T  of  unknown  magnitude  will  act  on  the 
underside  of  the  deckhouse,  and  a  similar  force  T  will  act  in  the  oppo- 
site direction  on  the  hull.   Equilibrium  requires  that: 

Nx  =  -  T,   Mx  -  -  aC^T  (la) 

N2  =  T,    M2  =  M  -  aa2T  •  (lb) 

Since  it  is  assumed  that  Navier's  hypothesis  is  valid  for  deckhouse 
and  hull  separately,  the  longitudinal  stress  distribution  across  a 
section  can  be  determined  for  the  deckhouse 

•  Nl   Ml 

Gi  =  A7  "  T2  xi  '  (2a) 

and  for  the  hull 

N2   M2 
°2   =  h.  '    h   ^  (2b) 

with  tension  stresses  counted  as  positive. 


At  the  junction  of  the  house  and  hull,  the  stresses  Q1  and  a_ 
must  be  equal,  with  x.  =  -  act  and  x„  ■  act, .   Using  Equations  1  and 
2,  the  value  T  is  found  to  be: 


T  = 


aVl 


M 


and 


I1I2   I1I2    2,  2   _     2  _  N 

——  +  -r —  +  a  (a_  I  +  a.  O 

A1     A„        2   112 


dT  a0t2Il 

dz   II    I  I2    2.  2  T     2  _  , 
~  +  "A—  +  a  (a2  ri  +  Ql  I2) 


(3) 


(A) 


T  was  defined  as  the  total  horizontal  shear  force  acting  between  the 

left  end  of  the  deckhouse  and  the  section  at  z.   The  total  shear  force 

in  the  structure  V  is  equal  to  dM/dz.   According  to  Equation  3,  the 

shear  T  at  the  ends  of  the  deckhouse  is  not  zero.   At  points  slightly 

away  from  the  deckhouse  ends,  there  is  no  deckhouse,  and  hence  T  =  0. 

This  means  that  in  addition  to  the  distributed  shear  force  dT/dz, there 

must  be  concentrated  shear  forces  T  and  T,  at  the  ends  of  the  deck- 

c      d 

house.   These  concentrated  shear  forces  cannot  exist  in  any  real  struc- 
ture.  Their  appearance  is  due  to  the  fact  that  shear-lag  effects  were 
neglected  when  Navier's  hypothesis  was  assumed  to  be  correct  for  the 
full  length  of  the  deckhouse.   These  forces  will,  in  reality,  distri- 
bute themselves  over  a  finite  distance,  presumably  over  a  distance  equal 
to  the  depth  of  the  deckhouse  as  depicted  in  figure  1.3. 


IIGURE  1.3 
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Vith  this  in  mind,  the  stresses  found  usinp  Equations  2  will  be 
incorrect  for  calculations  near  the  ends  of  the  deckhouse.   Accord- 
ing to  Saint  Venant's  theorem,  however,  the  effects  of  the  simplifica- 
tion will  not  affect  the  stresses  in  the  center  portions  of  the  struc- 
ture. 

In  simplifying  the  expressions  for  the  moments  and  direct  forces, 
the  following  notations  are  used: 

2   Al  A2 

I.  =  a"  ~~  <5> 

A      Ax  +  A2 

11  +  ^    ^ 
U  =  _--_  .  (6) 

2   a2  A 

The  term  p  is  a  measure  of  the  size  of  the  deckhouse  in  relation  to 

the  hull  and  is  referred  to  as  the  size  factor.   The  total  moment  of 

inertia  I  of  the  section  can  be  expressed  in  terms  of  I  ,  I-  and  by 

the  constant  I. 
A 

I  -  h   +  I2  +  IA  .  (7) 

With  the  use  of  notations  (6-7),  the  following  important  expressions 

are  derived: 

MIA    MIA     uCc^  -  ua2) 
N2  =  "  Ni  =  T  =  aT~  "  ~  (1  +  p)  (a2l^+  ya^)         (8a) 

MI 

M.  =  ~-  -  MI,  7r-T-w--T— I T-T  (8b) 

1  I     1  (1  +  y)(a2I1  +  pa1I2) 

MI  2 

2  I     2  (1  +  y)(a2Ix  +  ya^) 


These  expressions,  along  with  Equations  2,  could  be  used  to  determine 
a  and  a_  in  the  deckhouse  and  hull.   It  will  be  noticed,  however, 
that  Equations  8  contain  two  terms  the  first  term  being  the  value  of 
the  respective  N  or  M  if  the  conventional  bending  theory  were  applica- 
ble to  the  entire  section.   The  actual  stresses  CT  can  now  be  expressed 
as  the  sum  of  the  stresses  0  according  to  Javier's  hypothesis  and  a 
correction  (Aa) . 

a  -  cr>T  +  Aa  (9) 

with  the  Nav:  >:r's  stresses  (a,  J 

%--¥  '  .  ao) 

where  x  is  the  distance  from  the  centroid  of  the  entire  section, 

figure  1.1.   The  corrective  stresses  Aa..  and  Aa„  in  the  deckhouse 

and  hull,  respectively,  are 

AN    AM 
Aa.   =  ---  -  — -  X-  (11a) 

1    Al    h        1 

AN    A*<0 
Aa.  -  r--  _  *   y  (lib) 

z    A     I     I 

whore  AN  and  AN  are  the  corrective  portions  of  N  and  M  given  by  the 
second  terras  of  Equations  8. 

Although  the  model  in  this  analysis  was  simplified  in  assuming 
that  no  vertical  forces  act  between  hull  and  deckhouse  (i.e.,  zero 
deck  stiffness),  the  important  result  remains  that  the  stress  distribu- 
tion can  be  expressed  as  the  sum  of  the  stress  a..,  according  to  Navier's 

N 

theory  and  the  corrective  Aa.   A  typical  example  of  stress  distribution 
using  this  analysis  is  shown  in  figure  1.4. 
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General  Analysis  of  Two-r         ture 
In  this  analysis  the  assumption  is  made  that  any  relative  displace* 
tit  of  Che  deckhouse  with  respect  to  the  hull  is  resisted  by  the  inter- 
rl  vertical  forces  required  to  deflect  bulkheads  and  transverse  beams 
;  pporting  the  deckhouse.   External  vertical  loads  and  bouyancy  will 
(  ise  the  structure  to  deflect,  and  this  deflection  can  be  described 

the  displacements  y.  and  y?  of  the  center  lines  of  the  deckhouse  and 
t  LI,  respectively,  Mgure  1.5. 

DEFLECTED       CtMtlR  LIN  E  1 


A 


m—   C.G.  OF  0E.CKHOU5I  ^  \ 

~3- b-.-Jr^ivzPt 


-"" C.C.   OF   HULL 


HGURE  1.5 


Tl  stiffness  of  bulkheads  or  deck  beams  resisting  relative  vertical 
1  placements  is  assumed  to  be  constant  tor  the  full  length  of  the 
dfkhouse.and  the  magnitude  of  this  stiffness  is  given  by  a  spring 
-tant  K.   K  is  defined  by  Bleich  as  being  the  force  per  unit  length 


of  deckhouse  required  to  produce  a  relative  deflection  eoual  to  one 
unit  of  length.   The  vertical  reaction  between  hull  and  deckhouse 
will  be  therefore  equal  to  K(y  -  y  )  per  unit  of  length. 

Using  the  theorem  of  stationary  potential  energy,  the  differen- 
tial equations  for  the  deflections  y..  and  y„  can  be  obtained.   The 
potential  energy  U  consists  of  the  internal  strain  energy  and  the 
potential  energy  of  the  external  forces. 

U  "  2  '-1/2    [EI1  h    '  +  EV?.    "  +   EV«1V1     +  V2   ' 
+  K(yx  -   y2)2   -   2plVl  -    2p2y2]dz  +    [Vf'2V\j2 

-   [Sy2]*02  .  (12) 

U  will  be  a  minimum  if  the  variation 

6U  =  0  .  (13) 

Using  the  calculus  of  variations,  the  set  of  two  simultaneous  equa- 
tions are  derived. 

E(I1  +  VA)ylV  +  Kyl  +  aia2  EIAy2V  "  Ky2  =  Pl  (Ua) 

V2  EIAylV  "  Kyl  +  K(I2  +  ^A^T   +  Ky2  =  P2  (14b) 

p1  and  p„  are  the  distributed  loads  acting  on  the  deckhouse  and  hull, 
respectively. 

The  process  of  the  calculus  of  variations  also  furnish  the  boundary 
conditions  required  to  determine  eight  arbitrary  constants  which  will 
appear  in  the  general  solutions  of  the  differential  equations.   At 
z  =  +  1/2   and  z  =  -  1/2 
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y2   =   0  (15a) 

EClj  +  a*    lA)y1     +  01^2   EIA  y2"   =  0  (15b) 

aia2    LIAVl'   +   E(I2   +  a2    1t\)v2      =  "  M  (15c) 

E(I.    +  a2;    I  )y      '   +  a.  a.   EI  v"     =  0  (15d) 

1  1       A       1  1    Z  A    £ 

Hie  second  and  third  boundary  conditions  are  equivalent  to  M.  =  -  aa,T 
and  H_  =  M  -  acuT  as  presented  in  Equations  (1).   The  fourth  boundary 
condition  indicates  that  the  shear  forces  at  the  end  of  the  deckhouse 
must  vanish. 
3,   Solution  of  Differential  Equations  for  Constant  Voreent 

In  reference  to  the  free-body  diagram  of  tisure  1.6,  this  analysis 
neglects  the  loads  p..,  p«  and  the  shears  S  and  S  .   The  onlv  load  on 
the  model  is  M  =  M  ■  M  (i.e.,  constant  moment).   The  differential 
Equations  14  are  then  homogeneous,  and  since  the  problem  considered 
is  symmetrical  with  respect  to  the  origin  of  the  coordinate  z  axis,  the 

general  solution  contains  four  arbitrary  constants. 

2 
y  =  C  +  C?z  +  C-  sinyz  sinhyz  +  C,  cosyz  coshyz  (16a) 

2 
y2  =  C.  +  CLz   -  pC-  sinyz  sinhyz  -  UC,  cosyz  coshyz        (16b) 

where 

.4 


y  = 


4E  a2I1  +  ya1I2  K     J 


These  solutions  in  conjunction  with  the  boundary  conditions  furnish  the 
values  of  the  constants. 
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C2  -  "  2E7  (18*> 

C3  =  ^T  0^irrE("a2T~T17a1"lJ)'  (18b) 

x         U'^1  M 
CA  =  ^2"  a"+~vT)"Tra2T1"-f  "ua'ip"  *  (18c) 

C.  is  not  included  since  it  only  describes  rigid  body  motion.   The 

values  of  0..  and  ty.    are 

.    sin  u  cosh  u  +  cos  u  siiih  u  /nrk  N 

<j>  =  -. — . . (19a) 

1   sin  u  cos  u  +  sin  u  cosh  u 

i    cos  u  sinh  u  -  sin  u  cosh  u  mqk^ 

1   sin  u  cos  u  +  sinh  u  cosh  u 

u  »  li  »  S  \4fL.   _JL±JL_  (20) 

2         2     \|4E     a2Ix  +  ua^        *  UU; 

Introducing  the  equations  describing  the  deflections  y1  and  y^  into 
the  following  expressions  for  the  moments  and  direct  forces 

Mx  -  -  EI1  f1  (21a) 

M2  =  -  EI2  y2'  (21b) 

EIA 
Nx-  -  ::2  =  —  (al7l  +a2  y2  >,  (21c) 

the  following  values  were  found  at  the  amidhsips  section,  z  =  0. 

I.  I,         y(a     -  ua?) 

1  al       91    •  a        (1  +  y)(a2I1  +  uaxI2)  U  a) 

I 

M         =  >t  -i  _   o     MI-      7T-r—  w~T— r fT  (22b> 

1  I  1       1     (1  +  u)(a2I1  +  pa1I2) 

\                                          v2 
M„      =  M  ~  +  $.  mi.     7T-r— rr-~T~ r rr  •  (22c) 

2  I  1       2      (1  +  u)(a2   I.   +  ua1I2) 
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A  remarkable  result  is  that  If  a  comparison  is  now  made  between  Equations 
8  obtained  in  Section  1  and  Equations  22,  the  only  difference  is  that 
the  factor  $.  is  nov;  introduced  into  the  second  term.   Therefore,  the 
refined  theory  used  in  this  section  furnishes  a  similar  result  as  that 
obtained  in  Section  1  (zero  deck  stiffness)  but  modified  only  by  the 
multiplication  of  the  corrective  stress  Aa  by  a  factor  <J> 


o  =  o.,  +   <*>  Aa 

h  1 


(23) 


This  result  is  surprisingly  simple  in  that  the  deviation  of  the  stress 
distribution  from  Javier's  calculations  is  determined  by  a  non-dimen- 
sional factor  $- .   The  corrective  forces  AN.J ,  AN„  and  moments  AM-,  AM« 
are  computed  from  the  equations  given  in  Section  1 


|- i/JL 


-*•   I  *- 

I 


l/A 

—  2 


»c 


TTTTTTTTTT 


M 


.OAD5    P,  -*■  Pa 


HGURE  1.6 


^ '      Solution  of  Differential  Equations  for  Equally  Distributed  Loads 

In  this  final  analysis  Bleich  considered  only  the  case  of  equally 
distributed  loads  p.  and  p9  acting  on  the  deckhouse  and  hull,  respec- 
tively, where  the  moments  at  the  ends  of  the  deckhouse  are  set  equal  to 


.oo_ 


zero  (i.e..  M  =  M_  =  0).   Equilibrium  requires  that  the  external  shear 
forces  S  =  -S_  =  II '«  (p,  +  p9) .   The  moment  at  the  center  of  the  sec- 
tion  (1'ipure  1.6)  duo  to  the  loading  Is 

(P,  +  Po)   2 

Mp  -  -  ■■  3   -■■  I        .  (24) 

The  general  symmetrical  solution  to  Equations  14  are 

2 
y  =  C  +  C_z  +  C_  sin  yz  sinh  yz  4  C,  cos  yz  cosh  yz 

(pl    +  P2}      A  Pl 

4  ~Tf~z  +  rr+r:  (25a) 

2 
y_  =  C1  +  C?z  -  UC„  sin  yz  sinh  yz  -  yC^  cos  yz  cosh  yz 

(Px  +  P2)   4       P2 

+  -2— i- z  +  rr+TDic  •  (25b) 

The  boundary  conditions  expressed  in  Equations  (15)  furnish  the  values 
of  the  constants 

C2  "  "  jfe  <26a> 

1  U  $2  Mn 

C3  =  ^2"     (T^TTe^I"  +  pojlj)  (26b) 

C4  =  2y2     Tl  +  P)  E("a~I~  +  iTaii2)   •  (26c) 

and  where 

*    2       sin  u  sinh  u  /-,-.% 

j  =  —  — — .  (27a) 

2  .  u  sin  u  cos  u  +  sinh  u  cosh  u 

.i.    2       cos  u  cosh  u  ,0,,  N 

2   u  sin  u  cos  u  +  sinh  u  cosh  u 

The  values  of  the  constants  are  similar  to  those  in  the  preceding 

section,  and  the  remaining  solution  for  the  stress  distribution  at 
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z  =  0  is  the  same,  the  onlv  difference  being  that  a  no"  deviation 

factor  <t>0  replaces  $,  and  M  replaces  M. 
2  1      p 


N     2 


(28) 


Curves  of  $.  and  $9  are  plotted  as  a  function  of  u  in  figure  1.7 


i.o 


$ 


»  3 

ununF  1.7 


Since  u  is  a  function  of  K«  it  will  be  noticed  that  for  the  case  of 
zero  deck  stiffness  (K  =  0)  that  ?«  and  &«  are  both  eoual  to  1.0,  and 
bence  the  results  of  section  one  are  ^"f^  rmi  ^  T  =  Cw  +  A(7.   t,Then  K 
or  th°  dimensions  of  the  deckhouse  became  lar^e  and  bence  U  >  2,  the 


factors  '   and  5„  became  very  s   11  indicatin   that  "  •"'  ir's  i  c 
distributi    i   ipplicable  at  t:  .-   amidships  section. 
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cil  ::  -  EXPA.rsion  n)  'r  "r r 

Although  Bleigh's  analysis  dealt  primarily  with  the  stress  dis- 
tribution located  at  the  center  of  the  deckhouse  (z  =  0) ,  it  is  the 
intention  of  the  present  analysis  to  also  investigate  the  stress 
distributions  away  fror  the  center.  Although  Eleich  mentions  that 
the  relationship  for  stresses  along  the  deckhouse  will  be  in  the  form 
of 

o  =  a..  +   $(z)  La     ,  (1) 

no  expressions  for  <5  as  a  function  of  z  vere  provided.   Tbey  can, 
however,  be  readily  formulated. 
1.   Constant  llor.cnt  Loading 

lor  the  case  of  constant  moment  loading, the  expressions  for  the 

deflections  are: 

2 
y  =  C  +  C0z  +  C_  sinyz  sinhyz  +  C,  cosyz  cosbyz  (2a) 

2 
y0  =  C  +  C^z  -  yC~  sinyz  sinhyz  -  uC,  cosyz  coshyz  .      02b) 

The  second  derivatives  of  y1 ,  v?  T,ith  respect  to  z  can  be  introduced 

into  the  expressions  for  the  moments  M_  ,  Vf0,  and  N- ,  N„ . 

Mx  -  -  KIj  t'i     •  M2  -  -  n,  ■•;'  (3a) 

V  T 

Nl  =  -  N2  =  T"   (aiyl'  +  a2y2?)  (3b) 

yielding  the  follov/ing  results: 

"  MT 
1 
M.  =  -= [0.  cosyz  coshyz  -  ij/.  sinyz  sinhyz] 

MI1  Cl  +  u)  (a^  +  vafj  (Aa) 


!'     =  — =—  +   [$.    cosyz   c  -  ik    sinyz 


n 


— - 
fl 


•  p)(a2Ij  +  yc^Tj 


4\ 
HL   =  -  N.    =   ~ —    [  !-     cosyz   coshyz  -  \\t.    sir.yz  sinhyz] 

MI  y(a]   -  ua0) 

These  expressions  arc  exactly  in  the  fori  oi      Lons  22  of  Chapter  1 
except  that  a  net:  deviation. factor  now  e::.'  '    ich  that 

a   =  o\,  +  *x(z)  Aa  (5) 

where 

0,  (z)  =  [0.  cosyz  coshyz  -  iy.  sinyz  sinhyz]    .  (6) 

In  order  to  modify  the  argument  of  the  trigonometric  and  hyperbolic 

yic  z 

functions,  it  will  be  remembered  that  u  ■»  ±~  ,  and  hence  yz  =  2u  -  , 

*Tow  *  (z)  can  be  expressed  in  terms  of  a  non-dimensional  distance  from 

the  center  of  the  deckhouse. 

z         z  z         z 

$.  (z)  =  1'  cos  2u  r-  cosh  2u  r-  i|/-  sin  2u  r-   sinh  2u  t-  .     (7) 

With  the  aid  of  tORTILAU  comput  :r   tl  w.is  to  faci]  Li       Lculal  '^ns, 
values  of  $-  (z)  were  calculated  as  a  function  of  u  and  (z/ Z) .      A  family 
of  curves  for  different  values  of  u  is  presented  in  tigure  2.1.   In 
finding  the  bending  stress  distribution  away  from  the  centerline,  no 
additional  calculations  need  be  made  other  than  selecting  the  appro- 
priate deviation  factor  from  the  provided  curves  and  multiplying  it  by 
the  corrective  stress  Aa. 


o 


'-n 

_o 
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2.   Equal  1 v  Distributed  Loads 

tor  the  case  of  equally  distributed  loads,  the  equations  for  the 

deflections  are: 

2 
y  =  C  +  C_z  4  C,  sinyz  +  C,  cosyz  coshyz  + 

(pl  +  P2)z4      pl  (c   v 

4  T4— *  (1  +  y)K  C8a) 

2 

y.  =  C  +  C0z"~  -  pC-  sinyz  sinhyz  -  \iC,    cosyz  coshyz 

(p,  +  P0)z       UP, 

tollowing  the  same  procedure  used  in  the  preceding  section,  the  expressions 
for  the  moments  and  direct  forces  are: 

M  =  -H-A  -  [$  '  cosyz  coshyz  -  |/-  sinyz  sinhyz] 

u       •     (pl  +  P25   Tl   2 

M  T   H. r-  -  - —   --   z  CQt) 

pi  (1  +  vOCa^  +  ua^IJT      2      I  ua; 


M  I 


o 


M„  =  _L_±.  +  [$_  cosyz  coshyz  -  (Jl-  sinyz  sinhyz] 

p2  (?!  +  D2)      JJ      2 

MpX2    U^TOTa^Tl^l")"  2  I"  Z  (9b) 

M  IA 

:*     =  -   N,    -     °-        -    [\    cosyz    coshyz  -  ,'j0   sinyz   sinhyz] 
/  1         al  z  z 

MpTA     _^(CC1  "_^2) ^  Pji  I*     2 

•        a  (1  +  u)(ct2I1  +  ya1I2)     "  2a  I     2      '  Uc; 

At  first  glance  these  equations  do  not  appear  to  be  in  the  standard 
form  due  to  the  presence  of  an  additional  third  term.   The  first  and 
third  terms  of  each  expression  can  be  combined,  however,  such  that 


- 


N 


fei 

< 

Ll. 

Z 

o 

f- 


_ 


Vi    .(pi^p2)    h    2     Vl 

i -2 —  r  7  =  —  (10a) 

Hh        <P     +  p2)    I2     2       MZI2 

i 2 —  r z  =  t-  (10b) 

M   I,         (p.    +  p„)    I,      9        !1I 

.£  A  _        1    _„2_     A   ^2   =      z   a  (10c) 

al  2a  I  al 

where  M  is  the  moment  due  to  the  distributed  loads  located  at  a  dis- 
z 

tance  z  from  the  ccntcrline.   Equations  0  can  then  be  used  to  express 
the  stress  distribution  in  the  form  of 

a  =  a:;  +  ■;  (z)  Aa  (11) 

where 

z         z  z         z 

$_(z)  =  $-  cos  2u  t*  cosh  2u  » —  lp-  sin  2u  j   sinh  2u  r-  .     (12) 

It  must  be  remembered,  however,  that  the  Javier  stress  term  is  deter- 
mined through  the  use  of  the  moment  (M  )  at  the  location  In  question 
whereas  the  corrective  stress  term  (Aa)  employs  the  use  of  the  moment 

at  the  center  of  the  structure  (M  ). 

P 

As  in  the  previous  section,  computer  calculations  were  used  to 
find  &«(z)  as  a  function  of  u  and  (z/£) •   A  family  of  curves  for  differ- 
ent values  oi  u  is  presented  in  lipure  2.2. 

tor  the  case  of  constant  moment  loading,  the  value  of  0  (z)  at 
the  ends  of  the  deckhouse  (z/£  =  .5)  is  equal  to  1  for  all  values  of 
u.   The  reverse  is  true  in  the  distributed  load  case.   $~(z)  is  eaual 
to  zero  at  the  ends  of  the  deckhouse  for  all  values  of  u.   This  is  as 
expected,  since  the  moments  at  ends  of  the  deckhouse  must  vanish  in 
the  distributed  load  model,  and,  therefore,  there  can  be  no  corrective 
stress  term  at  that  location. 


- 


3.  Linear  Sunemositi  n 

In  reference  to  obtaining  a  tote]  solution,  Bl  Lch  did  not  mon- 
m  or  su    t  the  application  ^f    linear  superposition  to  hir-  two 
generalized  solutions,   ^n  the  other  hand,  there  is  nothing  to  Indi- 
cate that  the  amplication  of  linear  superposition  ir;  not  valid.   The 
simplifying  method  of  breaking  the  free-body  diagram  cf  fcigure  1.6 
into  a  constant  moment  part  and  a  distrj   Led  1   ' :    art 

id  shear  fc    N 
it  actual  loading,  shear,  and  1       moment         for 

the  model  when  the  two  loading  conditions  arc  then  summed  together. 
With  the  assumption  that  this  reasoning  also  applies  to  the  respec- 
tive stress  distribution  solutions,  I    t  tal  stress  distribution 
(O^)   will  be  expressed  as  the  linear  sum  of  the  stress  due  to  a 

constant  no:.:ent  (Cw)  and  the  stress  due  to  a  distributed  loading  (o  )  . 

M  ;: 

aT  =o^+o  (13) 

re 
aM  =  a,.  +  $.(z)AaM  (14) 

a  =  a.,  -i  $0(z)Aa   .  (15) 

p     ..     2 

Care,  however,  must  be  taken  such  that  the  appropripte  moments  are 

used  in  these  equations.  La      for  the  deckhouse,  as  an  example,  employs 

the  use  of  the  moment  y\ 

P 

MI.     u(a  -  \Kt0) 


A1  a   (1  4  ii)(a2I1  +  p^fT    P  (1  +  u)(a2I1  +  ^a1I2)  l 

as  described  by  Equations  11  in  Chapter  1.   Aa  is  different  only  in 

that  it  makes  use  of  the  constant  moment  M  in  lieu  of  !f  .   With  this 

p 

in  mind.  Equations  14  and  15  may  be  combined  such  that: 


L- 


M  x 

0 


"L 


=   "  i~  "  T~  +  :'i(z)  :":'  +  :2(z)  Ao  (16) 


(::  +  ::  )x 


aT  =  •  _  -   +  #.  (z)  AaM  +  $,(z)  Aa  (17) 

T        I  i       i     2      p 

In  later  chapters  this  assumption  of  linearity  v:ill  be  checked  by  com- 
parison of  results  on  similar  models  using  finite  element  methods  and 
Bleich's  method. 


CHAPTER  III  -  ICES  STPxUDL  METHODS 

The  Structural  Design  Language  (STPUDL)  developed  at  the  M.I.T. 
Civil  Engineering  Systems  Laboratory  is  a  series  of  computer  programs 
for  solving  structural  engineering  problems.   It  is  a  subsystem  of 
the  Integrated  Civil  Engineering  System  (ICES)  and  can  be  applied  in 
many  ways  to  a  wide  class  of  practical  problems.   An  operational 
system  has  been  implemented  on  the  IBM  system  /370.   '.-.Tien  entered 
into  the  computer,  the  STRUDL  language  Is  translated  by  the  ICES 
Command  Interpreter,  which,  for  each  command,  calls  upon  the  appro- 
priate STRUDL  program.   In  essence  it  is  a  language  with  which  an 
engineer  can  describe  a  problem,  its  solution  procedures,  and  ask  for 
results. 

Of  the  various  analyses  offered  in  STRUDL,  the  stiffness  analysis 
was  employed  for  the  present  investigation.  The  stiffness  analysis  is 
a  linear,  elastic,  static,  small  displacement  analysis. 

The  basic  principle  behind  finite  element  techniques  is  the 
replacement  of  the  actual  physical  problem  by  a  model  composed  of  a 
finite  number  of  discrete  elements  which  are  connected  at  their  nodal 
points.   Many  elements  are  available  which  have  different  geometric 
shapes  and  different  applications  for  the  solution  of  different  types 
of  problems,   tor  the  present  application,  the  primary  element  type 
used  is  the  'PSR'  element.   The  'PSR'  is  a  rectangular  element  with 
four  nodal  points,  one  located  at  each  of  the  corner  joints.   Used 
in  conjunction  with  a  plane  stress  problem,  the  nodal  unknowns  or 
displacements  are  U..  and  U9. 


U  =  a.  +  a9x  +  ay  +  a,xy 


2    5    6     7     8 
No  element  loads  are  available,  and  hence  all  loads  must  be  in-plane 
joint  loads.   Uhile  resulting  displacements  are  given  at  the  nodal 
points,  stresses  and  strains  are  given  only  at  tho  baricenter  of  each 
element,  tigure  3.1. 
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IGURE   3.1 

*  X 


The  material  specified  must  be  isotropic,  and  the  only  relevcnt  ele- 
ment constants  that  need  be  declared  are  the  foung's  modulus  and  the 
Poisson  coefficient. 

Another  element  type  used  was  the  "BPR'  clement.   It  is  usee! 
in  plate-bending  type  problems  and  is  similar  in  geometry  to  that  of 
the  'PSR'  element.   The  'CPU'  element  allows  distributed  type  loads 
in  the  out-of-plane  direction.   In  later  analyses  the  'BPR'  elements 
were  superimposed  on  the  'PSR'  elements  in  order  to  obtain  a  more 
refined  analysis. 
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tor  this  analysis  a  model  was  selected  vith  arbitrary  dimensions 
as  shown  in  Mgure  4.1.   Bulkheads  are  placed  at  equally  spaced  dis- 
tances of  20  feet  in  the  hull  section.   The  thickness  of  the  hull, 
box  girder  plating  is  .5  inch,  and  the  thickness  of  the  deckhouse 
plating  and  bulkheads  is  .25  Inch.   The  material  constants  inclu 
a  foung's  modulus  of  30  x  10   psi  and  a  Poisson  coefficient  of  0.3. 

The  mode]  I         15-foot  Iraft  wit!  a  correspon  - 

ing  hydrostatic  upward  distributed  force  of  17.143  tons/ft.   The 
internal  loading  is  arranged  so  as  to  provide  equi     urn  and  a 
resultant  symmetrical  loading,  figure  4.2.   Shear  and  moment  din- 
for  the  total  model  are  also  provided. 
!•   STRuT)L  Model 

Due  to  the  symmetry  of  the  loading,  it  is  necessary  onlv  to 
look  at  one-quarter  of  the  deckhouse  and  approprd  t   hull  structure. 
This  also  aids  in  reducing  computer  time  which  can  become  excessively 
In'gh  when  dealing  with  finite  element  methods.   The  planes  of  symmetry 
for  the  p  rob  1    r  i  tl   ~  ~  '        1      It!      '  r  lir    "    of 
the  model  in  1 igure  4.1.   The  one-quarter  structure,  broken  into  122 
'  PSP. '  elements,  is  shown  in  figures  4.3  to  4.5. 

Since  the  'PSR'  elements  only  alio"  in-plane  joint  loads,  it  was 
necessary  to  model  the  shear,  moment,  and  distributive  loads  into 
forces  acting  on  joints  or  nodal  points.   The  pressure  loading  acting 
on  the  bottom  of  the  hull  structure  (319.984  psi)  was  translated  into 
forces  acting  at  bulkhead  locations.   The  moment  at  the  end  of  the 
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HGURE  4.2 
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deckhouse  C'n)  ,;'ns  translated  into  forces  acting  on  the  end  joints 
of  the  hull  girder.   Using  the  notation  in  faction  3.3  of  reference 
6,  the  vertical  shear  flow  (")  at  the  end  of  the  structure  was  found 
using 

3  o   2Af  +  *[1  -  e")2] 
2  4e    3Af  4  a 

where 

o  -  sD   . 

The  shear  flow  ^m^v.v.q'1   a  parabolic  distribution,  tl     .inum  occurring 
at  the  neutral  axis  of  the  hull.   The  shear  flow  '.'as  then  translated 
into  vertical  forces  acting  at  nodal  points  at  the  end  of  the  hull  side. 
The  final  loading  is  depicted  in  Mgure  4.6. 

A  sanple  computer  program  for  this  loading  is  provided  in  Appendix  A. 
2  •   STRUDL  Results_ 

In  order  to  find  the  magnitude  of  any  distortion  due  to  the  applica- 
tion of  point  loads  on  the  structure,  the  variation  of  longitudinal 
stress  along  the  deckhouse  length  was  plotted,  Mgure  A. 7.   trom  these 
results  it  appears  that  some  distortion  occurs  near  the  deckhouse  end 
as  evidenced  by  the  abrupt  chanee  in  the  stress  pattern.   These  effects, 
however,  appear  to  be  concentrated  only  within  a  distance  of  10  feet 
of  the  deckhouse  end. 

As  previously  mentioned,  the  stress  output  in  STRUDL  is  given  at 
the  baricenter  of  each  element.   It  will  be  noticed,  however,  that  the 
baricenters  of  the  elements  adjacent  to  the  amidships'  location  lie  on 
points  5  feet  away  from  amidshiDS.   There  is  little  or  no  stress 
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variation  in  this  area,  and  hence  negligible  error  is  introduced  if 
the  stress  at  the  adjacent  elements  is  used  for  th   midships'  results. 

The  longitudinal  stress  distributions  at  ar .       and  at  a 
point  35  feet  from  amidships  are  shown  in  (figures  4.  i  and  A. 9,  res- 
pectively.  Pronounced  shear  lag  effects  are  evident  in  both  cases. 
This  is  primarily  due  to  the  relatively  low  deckhouse  length  (1)  to 
hull  breadth  (R)  ratio.   tor  this  particular  model  1/B  =  2.5. 

An  unexpected  result  is  the  discontinuity  in  the  loneitudinal 
stress  distribution  across  the  main  deck  occurring  at  the  junction 
of  hull  deck  and  deckhouse  side.   An  apparent  cause  for  this  dis- 
continuity is  the  nature  of  the  relative  deflection  of  the  hull- 
deckhouse  connection  with  respect  to  the  bull  side,  HgurcrA.10  and 
4,11.   The  longitudinal  stress  due  to  bending  is  proportional  to  the 

curvature  of  the  beam  (-) 

P 

a  =  -  E(i)  y  (2) 

where  p  =  radius  of  curvature. 

It  is  apparent  that  the  curves  of  deflection  of  the  main  deck 
at  side  and  the  hull-deckhouse  connection  as  presented  in  Mgure  A. 12 
are  quite  different  from  one  another.   Obviously  the  curvatures  of 
the  two  deflections  are  of  different  sign  or  magnitude  at  respective 
locations  and  hence  effect  the  stress  distribution  across  the  main 
deck  so  as  to  produce  the  results  shown  in  Mgurcs  A. 8  and  A. 9.   In 
order  to  support  this  reasoninp,  a  separate  STPUDL  program  was  run 
in  which  the  amidships' plane  was  held  rigid.   Although  this  was  an 
unnatural  condition,  it  had  the  effect  of  increasing  the  deck  stiffness 
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in  the     liate  vicinity  of  tl     Ldships'  plane.   The  hull  and 
deckhouse,  therefore,  tended  to  deflect  in  the  sane  manner,  and  the 
discontinuitv  was  virtually  eliminated  as  evidenced  in  Hgure  4,13. 

Another  area  of  interest  is  the  distribution  of  the-  shear  stress 
along  the  hull-deckhouse  connection.   ^uo  to  the  nature  of  the  ele- 
ments used,  it  was  only  possihle  to  obtain  results  of  the  shear- 
force  distribution  2,5  feet  above  tl   hull-deckhou   connectie^. 

However,  the  irregu]  r  |  if  tl        stress  as  shotm  in  Hg- 

ure  4.14  is  similar  Lr  ults  obtained  by  exoeriments  en 

physical  models  as  described  in  r^fer.^nco  2, 
3.  Bleich  Model 

The  initial  step  in  usinj  Bleicl  '  |  ig  to  bre  '  f",,r>  froe- 

body  diagram  of  fcignre  1,6  into  a  constant  moment  part  and  a  dis- 
tribu  '  .    Ltl;  refer  . 2 ,  it 

een  that  M  -  M-  ■  !*  ■  7,1^     I     ,  p.  +  p,  =  5.7 

ton    ' r '.  and   tha      S  ,    -■  =  285.7  tons,      "  '•'     tl  •  liter 

C    -v 

cross   tr  '    '    !      .  ]  ,  : I 

Le  pr    lure  to  u        Lu tion  i  -      i 

and  2  in   .  ix   to  arrive  at  a  solution. 

In  a  model  with  bulkheads,  however,  the  primary  difficulty 
is  involved  with  the  determination  of  the  deck  stiffness  or  spring 
constant  (K) .   Since  K  was  defined  as  the  force  per  unit  length  of 
deckhouse  required  to  produce  a  relative  deflection  equal  to  one 
unit  of  length,  it  is  apparent  that  the  value  of  K  will,  in  reality, 
vary  along  the  deckhouse  length  due  to  the  emplacement  of  structural 

bulkheads.   In  order  to  simplify  the  use  of  Bleich' s  method,  however. 
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an  average  value  of  K  i        leternd   I.  To  a       this  end  a 
symmetrical  portion  of  the  hull  structure  was  modelled  to  include 
a  bulkhead  and  attached  deck  and  hotton  plating,  figure  4.15.   This 
I-beam  type  structure  is  simply  supported  on  Its  ends  and  alio 
to  deflect  under  vertical  forces  applied  to  the  hull-deckhouse  con- 
nections.  In  the  analytical  procedure  used  for  determining  the 
forces  needed  to  deflect  the  hull-deckhouse  connections  1  inch,  it 
must  be  kept  in  r;ind  that  althougl    •  deflection  duo  to  shear  forcer, 
arc  negligible  in  most  cases,  in  short  deen  metal  beams  the  deflec- 
tion caused  by  shear  may  become  a  significant  portion  of  the  total 
deflection.   In  this  case  shear  contributes  a  major  portion  to  the 
total  deflection. 

The  equation  for  the  deflection  due  to  the  internal  moments  (•?„) 
was  calculated  through  the  use  of  singularity  functions.   The  expression 

for  the  deflection  due  to  shear  (jT  ^  was  obtained  through  the  use  of 

T 

the  method  of  unit  loads  as  desc:-'         :  :rence  7.   The  total 

deflection  (f_)  was  then  expressed  as  the  sum  of  i\.  and  f  and  set 

T  T 

equal  to  1  inch.   The  equation  was  then  solved  for  P  (the  distributed 
load  acting  on  the  hull-deckhouse  connections) .   The  expression  for 
K  is  twice  the  value  of  P 

K  =  2P  .  (3) 

4 
This  analytic  approach  to  K  yielded  a  value  of  1.86  x  10  psi  for 

the  present  model  under  consideration.   Sample  calculations  are 

provided  in  Appendix  B. 

As  a  further  check  on  K  a  STRUDL  program  using  'PSR'  elements 

on  the  model  presented  in  ligure  4.15  was  run.   Arbitrary  forces  (1) 
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were  applied  at  the  locations  Indicated,  and  the  resulting  deflec- 
tion at  the  hull-deckhouse  connection  noted.   The  force  was  then 
scaled  for  a  deflection  of  1  inch.   P  was  obtained  by  dividing  1 

by  the  width  of  the  flange  (240  inches).   Using  Equation  3.1'  was 

4 
found  to  have  a  value  of  1.53  x  10   psi.   The  disparity  between  the 

STRUDL  K  and  the  analytical  K  can  be  attributed  to  the  fact  that  in 

the  analytical  approach  the  deflection  calculations  apply  to  the 

neutral  axis  of  the  beam  only. 

Sample  calculations  for  the  solution  of  the  stress  distribu- 
tion at  amidships  using  Bleich's  method  is  provided  in  Appendix  C. 
4.   Comparison  of  Methods 

Comparison  of  Bleich  and  STRUDL  results  at  amidships,  15  feet 
from  amidships  and  35  feet  from  amidships,  are  shown  in  Mgures  4.16, 
4,17,  and  4.18,  respectively.   The  results  for  the  stress  distribu- 
tion in  the  vicinity  of  the  deckhouse  center  compares  more  favorably 
than  do  the  results  at  a  lc    Lon  35  £       y.  lu   to  ' 

fact  that  Bleich's  solution  neglects  shear-L:     .'Cts,  as  discussed 
in  Section  1  of  Chapter  1.   At  each  location   the  Bleich  stress  dis- 
tribution across  the  main  deck  and  hull  bottom  approximates  an  average 
value  of  the  STRUDL  distribution.   The  Bleich  distributions  along  the 
hull  side  produce  lower  peak  values  of  stress  as  compared  to  STRUDL, 
The  distributions  of  the  stress  along  the  deckhouse  side  and  top 
attributed  to  Bleich's  method  are,  in  general,  of  greater  magnitude 
than  that  of  the  STRUDL  results.   The  stress  distributions  for  both 
methods  were  transformed  into  moments  and  checked  against  equilibrium 
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conditions;  that  is,  the  values  of  the  moments  obtained  at  each  sec- 
tion corresponded  to  the  respective  moment  on  the  bending  moment 
diagram. 

Because  of  the  large  shear-lag  effects  found  in  this  particular 
model,  it  is  difficult  to  make  an  accurate  statement  as  to  the 
applicability  of  linear  superposition  to  Bleich's  methods  (which 
ignores  shear-lag).   It  appears,  however,  that  if  these  shear-lag 
effects  were  reduced, equilibrium  would  reouire  that  the  STPI'DL  results 
approach  Bleich's  solution.   This  will  he  later  checked  through  the 
use  of  a  model  whose  increased  1/B  ratio  significantly  reduces  the 
shear-lag. 
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CIIAPTLR  V  -  MODEL  WITH  INC   '   D  1/r 

In  this  analysis  the  sane  model  as  presented  in  Chapter  4  \;as 
used  except  that  all  dimensions  in  the  longitudinal  direction  <;cro 
increased  by  a  factor  of  fo  ^o  as  to  increase  the  1/B  ratio  from 
2.5  to  5.0.   In  order  to  r  '  '  Ln  a  b(  M    moment  Lmilar 

to  the  one  used  in  the  preceding  chanter,  the  distributive  loading 
was  reduced  by  -1.  factor  of  four:  i.e. ,  p.  +  p«  ■  1.42   tons/ft. 
fron  this  loading  S  =  ?Ti  ■  142.85,  md  M  ■  x'r   -  ML  ■  7.142.5  ft  I  ns. 
Due  to  the  increased  distance  bet'  ds,  the  value  of  the 

deck  stiffness  (".)  was  r     I  1  si. 

The  same  procedures  as  described  in  Chapter  4  "ere  used  in  the 
determination  of  the  STRUDL  and  Bleich  solutions. 
1.   STRUDL  Results 

The  longitudinal  stress  distributions  at  amidships  and  at  a  point 
70  feet  from  amidships  are  shown  in  Mgures  5.1  and  5.2,  respectively. 
The  most  notable  difference  in  these  results  from  that  of  the  preced- 
inp chapter  is  that  the  deckhouse  has  become  more  effective  in  the 
support  of  the  stress  distribution.   This  Is  ^.v.n   to  the  fact  that  tbe 
deckhouse  ler  tl  has  been  increased  to  20n  feet. 

The  discontinuity  in  the  stress  distribution  across  the  main 
deck  is  still  apparent  but  to  a  lesser  decree.   A  look  at  the  deflec- 
tion curves  of  Hgure  5.3  indicates  that  the  differences  in  curva- 
tures between  the  deckhouse  connection  and  main  deck  at  side  are 
not  as  great  as  in  the  preceding  model. 
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'  ant  rer.ult,  '       .  Ls  that      '  ear-lag  i  ffects 
have  been  reduced  to  a  negligible  anount. 
2.        rison  of  Methods 

Comparison  of  STRUDL  and  Bleich  results  at  amidships,  30  feet 
from  amidships  and  70  feet  from  amidships,  arc  shown  in  Mpurc  5. A, 
5.5,  and  5.6,  respectively.   Again,  the  difference  between  solutions 
grows  progressively  worse  as  one  travels  away  from  amidships  or  the 
center  of  the  deckhouse  structure. 

The  most  important  result  of  this  analysis       t  at  the  ami 
ships  section  (where  Bleich's  solution  is  most  applicable), the  two 
solutions  are  almost  identical.   The  argument  that  linear  superposi- 
tion is  applicable  to  Bleich's  two  reneral  solutions  is  therefore 
confirmed  by  these  results.   It  is  .     Ly  evident  that  if  the  shear- 
lag  effects  were  completely  done  away  with,  the  two  solutions  would 
correlate  exactly,  except  for  the  discontinuity  at  the  hull-deckhouse 
connection.   It  oust  be  kept        I,  1   ever,  that  I  '    ays  nothing 
of  the  applicability  of  Elcich's  total  solution  for  use  as  a  design 
tool.   In  real  life  structures  shear-lag  effects  are  inevitably 
present  to  some  degree  as  exemplified  by  all  the  STRUDL  results  thus 
far  presented,   further  discussions  on  this  topic  are  presented  in 
later  sections  after  all  model  results  have  been  presented. 
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CHAPTER  VI  -  MODEI  WITHOUT  BULKHEADS 

In  order  to  investigate  the  effects  on  the  stress  distribution 
due  to  a  more  flexible  deck,  the  bulkheads  of  the  model  presented 
in  Chapter  4  v/cre  eliminated.   .Ml  other  parameters  such  as  loading, 
lengths,  boundary  conditions,  etc.  remained  unchanged. 

tor  the  STRUDL  model,  'PSR1  and  ,BPP.'  elements  were  super- 
imposed upon  ti'     Le  structure-  so  as  to  allow  out-of-plane  loads. 
Due  to  the  nature  of  th<  Ltions,  "'  ich  allows 

a  considerable  degree  of  flexibility,  the  program  was  not  acceptable. 
It  is  also  noted  here  that  for  models  with  bulkheads  and  using  both 
'PSR'  and  ' BPR'  elements,  the  costs  were  prohibitive  due  to  the 
increase  in  the  number  of  elements.   In  order  to  obtain  reliable 
results  for  the  present  analysis,  only  'PSR'  elements  were  used.   The 
loads  on  the  hull  bottom  were  attached  to  transverse  stiffeners  in 
the  same  locations  where  the  bulkheads  had  once  been. 
1.   STRUDL  Results 

The  longitudinal  stress  distributions  at  amidships  and  at  a 
point  35  feet  from  amidships  ire  shown  in  Figures  6.1  and  6.2, 
respectively.   The  most  notable  result  is  in  the  degree  of  stress 
reversal  occurring  along  the  hull  and  deckhouse  sides.   Because  of 
the  relatively  flexible  deck,  the  deckhouse  has  become  less  effec- 
tive in  the  support  of  the  stress  distribution. 

Shear-lag  effects  are  as  pronounced  as  in  the  original  model; 
however,  the  discontinuities  at  the  hull-deckhouse  connections  are 
not  as  great.   Figure  6.3  supports  the  diminution  of  discontinuity 
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as  well  as  h      thting  the  increase  in  relative  vertical  displace- 
ments between  the  hull-deckhouse  ronnoction  and  main  deck  at  side. 

The  shear-stress  distribution  2.5  feet  above  the  hull-deck- 
house connection  is  shown  in  Bigure  6.  . 
2-   Blpich  Model 

Tho  absence  of  bulkheads  facilitates  the  analytical  approach 
to  the  determiri  Lion  of  K.   '  transverse  slice  of  tho  main  'Vr' 
can  be  assumed  to  lip  a  qirnl'1  slender  sunnorted  r>n  Its 

ends.   This  beam  is  subject  rn   vertical  VH.  rtinns  duo  to  distril  ted 
loads  (P)  so  is  to  deflect  tho  hull  deckhouse  connect ' 

pie  be;   ief]  etion  t1    •.     value  of  K  w   found  to 
1.736  psi.  r     LI  whei  for 

all  other  bul     d  models. 
3 .  Co  ''.'tl-.ods 


Comparison  of  Bleich  results  at  amidships,  15  foot 

from  amd        nd  35  feet  fr        hips,  ar        •   ■  •      C>  .5 , 
6.6,  and  6.7,  r<      Ively.  Again,  the  difl  Lutd 

grows  pr:  .    as  one  travels  t 

Ihe  Bleich  distributj      long  the  hull  side  produce  lover  peak  values 
of  stress  as  compared  to  STRUDL,  whereas  tho  distributions  of  stress 
along  the  deckhouse  top  and  side  are,  in  general,  of  greater  magni- 
tude than  that  of  the  STRUDL  results.   In  each  case  the  Bleich  stress 
distribution  across  the  main  deck  and  hull  bottom  approximates  an 
average  value  of  the  STPvUDL  distribution. 
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C01ICLUSI0NS  ANA  PvECODt]   '  CONS 

The  principle  of  lincnr  superposition  can  be  applied  to  Bleich*s 
two  general  solutions  in  order  to  obtain  the  total  stress  distribu- 
tion solution.   Because  Bleich's  method  neglects  shear-lap,  the 
use  oi  this  i  ethod  d  .\y   be  of  limited  significance  when  dealing  v;ith 
real-life  structures.   As  the  ]     tio  incrc     .  however,  Bleich's 
solution  becomes  r.ore  applicable  for  the  center  portions  of  the 
deckhouse-hull  structure . 

cause  of  the  relatively  favorable  conparison  of  Bleich  and 
STRUDL  results  in  the  center  portions  of  the  deckhouse-hull  area  for 
all  models  considered,  there  is  some  indication  that  Bleich  methods 
may  have  some  application  for  design  pruposes.   The  rather  simple 
analytical  methods  which  are  employed  for  the  determination  of  a 
stress  solution  is  also  a  factor  in  consideration.   All  of  the  results 
of  the  present  analysis  indicate  the  Bleich  method  is  conservative  in 
reference  to  the  stress  distribution  over  the  deckhouse  proper  and 
that  it  supplies  an  approximate  average  stress  across  the  main  deck. 
The  onlv  shortcoming  is  that  the  Bleich  solution  underestimates  the 
value  of  the  stress  at  the  junctions  of  the  hull  side  with  the  main 
deck  and  hull  bottom,   tor  example,  in  the  model  with  bulkheads  and 
1/B  =  2.5,  the  value  of  the  stress  at  the  main  deck  at  side  is  almost 
20  per  cent  less  than  the  STRUDL  result.   As  the  deckhouse  length  is 
increased  so  as  to  diminish  the  effects  of  the  shear-lag,  Bleich's 
solution  converges  towards  the  finite  element  solution,  as  evidenced 
in  fcigure  5.4.   This  supports  the  conclusion  that  Bleich's  method  is 
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morc  strongly  applicable  to  structures  having  relatively  long  deck- 
houses . 

In  order  to  confirm  the  finite  c.       results,  especially  in 
the  area  of  the  discontinuity  occurring  at  the  hull-deckhouse  connec- 
tion, it  is  recommended  that  a  physical  model  be  built  similar  to 
the  mathematical  models  presented  in  this  analysis.   Strain  and 
deflection  analyses  under  similar  loading  conditions  could  then  be 
used  for  verification,  and  augmentation  of  results  obtained  through 
STRUDL . 
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APPENDIX  A  -  SAMPLE  STRUDL  PR0G1 

The  following  is  a  STRUDL  program  for  finding  element  stresses 

and  nodal  displacements  for  a  model  with  bulkheads  spaced  20  feet 

apart . 

STRUDL  'TEST  1'  'HULL-DECKHOUSE  IIJTERACTU  ' 
$  ONE  OUARTER  STRUCT!" 
TfPE  PLANE  STRESS 

DEBUT,  ALL 
DUMP  1  ] 

joi::t  ar  \         .  ::s 

1  o.  o.  o.  s 

2  5.  o.  o.  s 

3  10.  0.  0.  S 

4  15.  0.  0.  S 

5  20.  0.  0.  S 

6  0.  10.  0.  S 

7  5.  10.  0. 

8  10.  10.  0. 

9  15.  10.  0. 

10  20.  10.  0. 

11  0.  20.  0.  S 

12  5.  20.  0. 

13  10.  20.  0. 

14  15.  20.  0. 

15  20.  ::.  . 

16  0.  30.  0.  S 

17  5.  30.  0. 

18  10.  30.  0 

19  15.  30.  0. 

20  20.  30.   . 

21  0.  40.  0.  S 

22  5.  40.  0. 
2  3  10.  40.  0. 

24  15.  40.  0. 

25  20.  40.  0. 

26  0.  50.  0.  S 

27  5.  50.  0. 

28  10.  50.  0. 

29  15.  50.  0. 

30  20.  50.  0. 

31  20.  0.  10.  S 

32  20.  10.  10. 

33  20.  20.  10. 
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34  20 

.  10 

35  20 

.  40.  10 

36  20 

,  50.  10. 

37  20 

.  0.  20.  S 

38  2; ) 

,  10.  20. 

39  20 

20.  2''. 

40  20 

30.  20. 

41  20 

40.  20. 

42  20. 

50.  20. 

43  0. 

0.  10.  s 

44  5. 

0.  10.  s 

45  10. 

0.  10.  s 

46  15. 

0.  10.  s 

47  0. 

20.  10.  S 

43  5. 

. 

49  10. 

20.  10. 

50  15. 

20.  10. 

51  0. 

40.  10.  S 

52  5. 

40.  10. 

53  10. 

40.  10. 

54  15. 

40.  10. 

55  0. 

0.  20.  S 

56  5. 

0.  20.  S 

57  10. 

0.  20.  S 

58  15. 

0.  20.  S 

59  0. 

20.  20.  S 

60  5. 

20.  20. 

61  10 

62  15. 

20.  20. 

63  0. 

40.  20.  S 

64  5. 

4  . 

65  10. 

40.  20. 

66  15. 

40.  20. 

67  20. 

0.  30.  S 

6S  20. 

10.  30. 

69 

in   "in 

70  20. 

30.  30. 

71  20. 

40.  30. 

72  20. 

50.  30. 

73  15. 

0.  30.  S 

74  15. 

10.  30. 

75  15. 

20.  30. 

76  15. 

30.  30. 

77  15. 

40.  30. 

78  15. 

50.  30. 

79  10. 

0.  30.  S 

80  10. 

10.  30. 

81  10. 

20.  30. 

82  10. 

30.  30. 

33  10. 

40.  30. 

84  10. 

50.  30. 

85  5. 

0.  30.  S 

36  5. 

10.  30. 
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37  . 

5. 

1. 

30. 

83  . 

5. 

30. 

30. 

89  . 

rj. 

40. 

30. 

90  . 

5. 

50. 

30. 

91  l 

3. 

0.  30.  S 

92  ( 

3. 

10. 

30.  S 

93  ( 

3. 

20. 

30.  S 

94  ( 

3. 

30. 

30.  S 

95  ( 

3. 

40. 

30.  S 

96  ( 

3. 

50. 

30.  S 

97  : 

). 

50. 

35.  S 

98  ! 

35. 

99  : 

LO. 

50. 

35. 

100 

10 

.  40 

.  35. 

101 

10 

.  30 

.  35. 

102 

10 

. 

.  35. 

]03 

10 

.  10 

.  35. 

104 

10 

.  o. 

35.  S 

105 

0. 

50. 

40.  S 

10G 

5. 

50. 

40. 

107 

10 

.  50 

.  40. 

108 

10 

.  40 

.  40. 

109 

10 

.  30 

.  40. 

110 

10 

.  20 

.  40. 

111 

10 

.  10 

.  40. 

112 

10 

.  0. 

40.  S 

113 

10 

.  50 

.  45. 

114 

10 

.  40 

.  45. 

115 

10 

.  30 

.  45. 

116 

10 

.  20 

.  45. 

117 

10 

.  • 

.   _>. 

118 

10 

.  o. 

45.  S 

119 

5. 

50. 

- 

120 

5. 

40. 

45. 

121 

5  . 

30. 

45. 

122 

5. 

20. 

123 

5. 

10. 

45. 

124 

5. 

0. 

45.  S 

125 

0. 

50. 

45.  S 

126 

0. 

40. 

45.  S 

127 

0. 

30. 

45.  S 

128 

0. 

20. 

45.  S 

129 

0. 

10. 

45.  S 

130  0.  0.  45.  S 

JOINT  RELEASES  tORCE  f  Z  MOMENT  X 

6  11  16  21  26  47  51  59  63  92  93  94  95  96  97  105  125  126  127  123  129 

JOINT  RELEASES   *0RCE  X  Xfo:XY,T    if 

2  3  4  5  31  37  44  45  46  56  57  58  67  73  79  35  104  112  118  124. 

JOINT  RELEASES  tORCE  Z 

2  3  4  5  31  37  43  44  45  46  55  56  57  53  67  73  79  85  91  104  112  118  124  130 
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33  TO   4^ 
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'  »E 
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.  TfPE 

Til]              -0.50 

72  TO   91 

.  TfPE    '???,'  , 

\;:ss  o.50 

92  TO   12 

!2   TfPE    'PSR' 

.25 

CONSTA] 

E  30000000.  ALL 
POISSON  0.3  ALL 
G  11538462.  / 
UNITS  KIPS 

DINC  'TOTAL' 
$  TOTAL  LOAD  INCLUDES  SHEAR,  MOMENT  AND  PRESSURE 
.in:     \DS 
V  PRESSURE  LOADS 
5  *ORCE  Z  -{  . 


2   3  4   11   15   21  25  *ORCE  Z  - 

12  13  14  22  23 

$   SHEAR  AND  !       ...    LOADS 

:  26. ( 

27   23   29  iORCE   f  53.3344 

30  10RCE    f  68.1072   Z   49.28 

36  iORCE   f   35.5936  Z  110.712 

42  IORCE   f  -35.5936   Z   110.712 

72  1-ORCE   f  -68.1072   Z   49. 

78  84   90  fcORCE   f  -53.3344 

96  tORCE   f  -26.6672 

UNITS  tEET  KIPS 

print  r 

STIllNESS  ANALfSIS 

UNITS   INCHES   POUNDS 

LIST  DISPLACEMENTS   STRESSES  Al 

A* 
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appendix  d  -  calci     ns  ion  determinik 

The  following  is  an  analytical  approach  to  the  determination 

of  the  spring  constant  (!'.)  for  a  model  with  bulkheads  spaced  20  feet 

apart . 

Basic  .'"•••  :\cl  at ure 

A    -  area  of  beam  cross  section 

K'    -  factor  depending  on  shape  of  bean  cross  section 

P    -distributed  load  at  hull-deckhouse  connection 

V    -   vertical  shear  due  to  actual  forces 

v    -  vertical  shear  due  to  load  of  one  pound  acting  at  the 

section  where  the  deflection  is  to  be  determined 

I..        -      deflection  due  to  internal  moments 
M 

f    -   deflection  due  to  shear 
T 

f_   -   total  deflection 

x    -   distance  along  length  of  beam 

<  >     indicates  singularity  functions 

reflection  Duo  to  M    rit  Onlv 
2 

df>-5  1  1 

EI  t-t  =  240[P  -   P  <x  -  120>*  -  p  <x  -  360>  ] 
dx z       x 

3  . 

EI  ?M  =  240P[  ~-  -  <x  -  120>  -  <x  -  360>  +  C,x  +  C, 
M         6       b  ~-6         1     2 


B.  C.  I     =  0  when   x  =  0,   480' 

E  =  30  x  106  psi 
I  =  8,748.005   inA 


7      =    -    (2.1   '"    :■:    Tf5)P      at    >:   =    120" 


Tinfi  action    Duo    to    '  ' 


1       f  Vv 
f      =    -    ,-77-  /    TTT     dx    (method   of    unit    loads) 
T  !.'         AC. 


T 


1  P(2j 

AG 


K,  m wcLjLXIQ  _+_^r 

6(2  4   12m  -f  ?.5r-2  +  15m  )  4  y(11  +  66m  4  135m2  4  90m3)  4 


v/here 


m 


30mn    (1  4  m)   4    '/.-••    (8  4  9n) 


2b  cf       b 

TT    n  =  h 


" 


I    BEAN    CROSS    SECTION       h 


-I 


e-t 


UJ 


t.  =  0.5"  ;  t  -  .25"  :  h  =  360",  b  =  240" 

f  M 


m  =2.66  :  n=  .666 
K'  =  .262 


It  is  customary  to  assume  that  or.lv  webs  of  structural  shapes  such 

as  channels  and  I  beams  resist  shearing  stresses  because  shear  stresses 

are  crall  in  flanges. 

t     =  -  (10.55  x  10"5)P 


1  T 


=   1 


»  » 


- 


1 " 

p 

K 
K 


=     -  P(10.76  x  10     ) 

=      .93  x   in'4 

=      2P 

L 

=      1  .86   x  10     -«?i 
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appt'^Uv  r  -   SA.VPLE  PJT'T""  CALCU1  MTONS 

Th^  following  is  a  solution  for  the  stress  distribution  at 
amidships  for  a  model  with  bulkheads  ^t^i  °0  frr>t  apart. 

B-"?ir      elatur 

a         -   -1 -■•  ^ •  - -, r- r>  v^t-uppr,  Cp^«-o-^^   f   >-      of  hull  and 

deckhouse  sect 
A-         -   cross  sec'  '    L 
*  ^  cr 

-  to  tal  '         31  t  i 

-  f;  T 
:  DTticnt  of  inertia  of  dec'  '    !  cross  section 
moment  of  j       of  hull  cross  section 
spring  const 
length  of  deckhouse 
constant  moment  part  of  loading 

-  moment  in  midship  section  due  to  the  loads  p.  J 
correction  of  direct  force  acting  on  deckhouse 
correction  of  direct  force  ac~     n  I  ull 
correction  of  moment  actin  on  '  "'  use 
correction  of  moment  acting  on  hull 
equally  distributed  loa 

u         -  parameter  described  by  dimensions  of  bull  and  deckhouse 

and  ot  stiffness  factor 
x-         -  vertical  distance  from  center  of  gravity  of  deckhouse 

cross  section 


I 

h 

h 

h 

". 

1 

:: 

M 

? 

-i 

■  -•  i 

1 

-": 

(p1  + 

P 

x„        -  vertical  distance  fron  center  of  gravity  of  hull 

cross  section 
x        -  vertical  distance  from  center  of  gravity  of  entire 

cross  section 
r*  -  ratio  of  the  distance  of  center  of  pravitv  of  dec' 

from  main  deck 
a.       -  ratio  of  the  distance  of  center  of  gravity  of  hul 1  from 

main  dec 
Y        -  parameter  described  by  dimensions  of  hull  and       use 

and  of  stiffness  factor  K 
\i  -      size  factor;  measure  of  the  size  of  the  deckhouse  in 

relation  to  hull 

CTW        -   stress  solution  for  constant  I 

.1 

cr        -  Navier*s  stress 
N 

0  -  stress  solution  for  distributed  load 

0"  total 

Act  -  corrctive  stress  for  dc   '   ~e 

AcJ.  -  corrective  stress  for  hull 

0  -  deviation  factor  used  in  constant  moment  solution 

#~  -  deviation  factor  u^ed  in  distributed  loading  solution 
Significant  Values  of  "athemntical  Model 


a  =  306"   A.  =  150  in2  A£  =  840  in2 

I  =  534,600  in4  T7  =  19.^40,000  in4 

M  =  M  =  19?.  x  106  in  lbs 

a  =  .4118  a2  =  .5882 

K  =  1.53  x  10A  nsi 


Bleich  Cm  cl^tionr 
"■      ^'1    A?  4 

I.   =     a'    ,        -v       =     11,917,309   in 

K  A1   +    V, 

T     =      I.   4  I,  +  I.   ■      V  ,891,909 
J  2  A 

•  Ti !  j  Ti 

I     +  a     1 
Because   in   this  model   x'  -  M  ,    (      )..     ~      ( \c) 


]  a       ?1  +  u)(a'   "i"  "h 


\N.      -      1«.9   x   10"    It    . 

1  1    (1    •!-    :  ;  (    .,7.    ■'  ..) 


2  2      (1  +  u)  (a„I.    - 


"'1  "1 


=12        J    •     .- 


1 


(dt_^        use  top)      =     2167  psi   ;      x-   -  5< 


Aa.    (main   deck)      =     -    344   psi    ;      x     =  -   126" 


AJ2  =     Af  "    if  X2 

Acf„  =     -    225   -    3.44   x0 

La2  (main  deck)    =     -    344  psi;      x2   =   ISO" 

Aa2  (hull  bottom)      =   394  psi;      x„  =  -   180" 


=  It  =  TL 


'  N 


■  i  ,. 


1  +  P 


S  '•:  oul. 


K"1T2 


u  =  1.785 


fror.  prap'as     =  .274,   >„  =  .  3C5 

("  l  M  ) 
c..    -    --   x 

OL5    (deckhouse    top)      =     -   3770 
(1  all  bottoi  )  =     27: 


x  = 


y 


313. 64" 

- 

-226.36" 


aT  =  a,.  +  a 
I         P 


oT  -  a.,  4  *  (z)  Cr),.  4  b  («)  (Act) 


at   z  =  0   ^(z)  h  ♦i  .   &2(z)  =  02 


linal  Solution 


aT  (deckhouse  top)  =  -  2390  psl 


," 


2147  .  ' 


a_  (1  nil  "  ■-■ 


-  L- 
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